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Summary 

'The  layered  system  is  composed  of  elastic,  viscous  and  acoustic  layers 
bounded  by  upper  and  lower  acoustic  half-spaces.  Uniform  subsonic  flows  may 
be  present  in  the  acoustic  layers  and  half-spaces.  The  method  of  dynamic 
stiffness  coupling  is  used  to  obtain  the  system  dynamic  stiffness  matrix  which 
connects  'spectral'  displacements  and  external  stresses  at  the  interfaces.  Par- 
field  sound  radiation  and  the  wavenumbers  uf  free-waves  are  found  from  this 
matrix  relation.  Numerical  results  include  wavenumber-frequency  plots  which 
exhibit  the  characteristics  of  instabilities. 
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Cartesian  coordinates 
spherical  coordinates 
etc.  displacements 

etc.  transforms  of  displacements 

Fourier  transform  wavenumbers 

Lam^  constants 

density 

pressure  in  fluid 
radian  frequency.  =27rf 
sound  velocity  in  fluid 
acoustic  wavenumber,  =2rrf/c 

wavenumber  of  longitudinal  waves  in  elastic  solid. 
=w/Cl=wV[(  X+2m)  /p] 

wavenumber  of  shear  waves  in  elastic  solid. 
=w/c-(-=w/  Vt  p/p) 

=  V(k2(i-ceM/k]2-a2-;8M  .  Im('y)>0 

=Vf k(_2-a2_^2) ,  im(-yL)>0 

=V(  k,|.2_flj2-^2)  im(7T)>0 

flow  velocity  in  fluid 
Mach  number.  =V//c 
thickness  of  layer 
=V[(x-xo)  2+(y-yo> 
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Dirac  delta  function 


INTRODUCTION 


Pestell  and  James  [1]  have  described  the  procedure  necessary  to  calculate 
the  far-field  sound  radiation  of  layered  media  composed  of  solids  and  fluids 
excited  by  time-harmonic  point  forces.  They  use  the  Lam^  decomposition  of  a 
vector,  viz. 

=  grad('P)+curl(j^)  (1.1) 

together  with  the  gauge  condition 


div(f)  =0  (1.2) 

to  faciiitate  solution  of  the  differential  equations  satisfied  by  the  displacements  of 
the  elastic  and  viscous-fluid  layers.  The  method  of  dynamic  stiffness  coupling  is 
then  used  to  assemble  a  system  matrix,  from  finite  element  matrices,  that  relates 
interface  'spectral'  stresses  and  displacements.  The  far-field  pressure  in  the 
upper  half-space  is  simply  related  to  the  normal  'spectral'  displacement  of  the 
uppermost  interface.  Fortran  programs,  written  in  single  precision  complex 
arithmetic,  have  been  used  successfully  to  calculate  far-field  sound  radiation  and 
free-wave  dispersion  plots. 

The  analysis  contained  herein  is  based  substantially  on  the  work  of  Pestell 
and  James,  but  important  differences  are  present.  First,  alternative 
decompositions  of  the  elastic  and  viscous  displacement  fields,  viz. 

u=grad(  F)  +  ( 3Q/3y-3 2h/3x3z.  -3Q/3x-3 2H/3y3z.  32H/3x2+32H/3y2)  (1.3) 

are  used  which  simplify  the  programming  and  effectively  halve  computation  times. 
Secondly,  because  of  a  recent  interest  in  the  effect  of  fluid  flow  on  wave 
propagation  in.  and  sound  radiation  by.  a  'thin'  plate  [2],  it  was  decided  to 
include  uniform  mean  subsonic  flows  in  the  acoustic  layers  and  half-spaces. 
Finally,  because  numerical  ill-conditioning  was  expected  at  the  low  frequencies  of 
interest  to  a  stability  analysis,  it  was  felt  necessary  to  develop  Fortran 
subroutines  to  simulate  double  precision  complex  arithmetic  on  a  PDP-11 /34A 
computer,  whose  arithmetic  word  length  is  32  bits. 

The  description  of  the  problem  is  given  in  Section  2  and  the  'spectral' 
dynamic  stiffness  matrices  of  the  elastic  and  viscous-fluid  layers  are  derived  in 
Sections  3  and  4.  The  'spectral'  dynamic  stiffness  matrix  of  an  acoustic  layer 
with  a  uniform  mean  subsonic  flow  is  given  in  Section  5.  and  the  acoustic  half¬ 
spaces  with  a  mean  flow  are  considered  in  Section  6.  The  Fortran  programs, 
particularly  the  subroutines  used  to  simulate  double  precision  complex  arithmetic, 
are  discussed  in  Section  7.  Finally,  some  numerical  results  are  discussed  in 
Section  8:  these  results  include  wavenumber-frequency  plots  of  layered  media 
with  flow,  which  are  expected  to  exhibit  the  characteristics  of  Instabilities. 


2.  PROBLEM  FORMULATION 
a.  General 

The  layered  system  is  composed  of  an  arbitrary  number  of  elastic,  viscous 
and  acoustic  layers  sandwiched  between  upper  and  lower  acoustic  half-spaces. 
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The  elastic,  viscous  and  acoustic  media  satisfy  respectively  the  exact  linearized 
equations  of  elasto-dynamics,  visco-dynamics  and  acoustics.  Uniform  one- 
-  dimensional  flows  may  be  present  in  the  acoustic  layers  and  half-spaces.  The 

1  layered  system  is  excited  either  by  time-harmonic  stresses  applied  at  the 

interfaces  or  by  a  time-harmonic  point  source  located  in  the  lower  half  space. 
The  time  factor  exp(-iwt)  is  omitted  from  all  equations  and  the  geometry  is 
shown  in  Figure  1 . 

I  b.  Fourier  Transforms 
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It  is  convenient  to  represent  the  variables  as  Fourier  transforms.  For 
example 

•(■CB 

Ux(x.y.z)  =  ( 1 /47r2)  j  J  UjjC  a. /3.  z)  expdax+ijSy)  dad>8  (2.1) 


where  the  'spectral'  displacement  u^  is  given  by  the  inverse  transform 

/3.Z)  =  j  j  Ux(x.  y.  z)  exp(-iax-i>8y)  dxdy  (2.2) 

-OB 

The  Fourier  transforms  facilitate  solution  of  the  differential  equations  by 
reducing  them  to  algebraic  equations. 

c.  Spectral  Dynamic  Stiffness 


The  'spectral'  dynamic  stiffnesses  of  the  layers  and  half-spaces  are 
matrices  which  relate  'spectral'  stresses  or  pressures  to  spectral'  displacements 
at  the  interfaces.  These  matrices. 

(S(a./3)1.  tSyCa.)?)].  .[F(a.iS)l.  (Fy(a./S)]  and  [FL(a./S)), 

6x6  6x6  2x2  1x1  1x1 

for  the  elastic,  viscous  and  acoustic  layers,  and  upper  and  lower  half-spaces, 
are  obtained  in  Sections  3-6. 

The  assembly  of  the  element  matrices  to  form  the  system  dynamic  stiffness 
matrix  relation 


lZ(a.j3)]lU{a.^)]  =  IE(a.>3)] 


(2.  3) 


is  a  standard  finite  element  procedure  that  reflects  continuity  of  displacement  and 
equilibrium  of  stresses  at  the  interfaces.  Computer  implementation  is 
straightforward.  The  matrix  (U(a.^)]  is  a  column  matrix  of  interface 
displacements  and  [E(a.)3))  is  a  column  vector  of  externally  applied  'spectral' 
stresses  which  are  defined  to  be  positive  when  acting  in  the  positive  direction  of 
the  coordinate  axes.  With  reference  to  Figure  2.  this  sign  convention  requires 
changes  of  sign,  before  assembly  of  the  elements,  to:  (i)  the  last  three  rows  of 
the  elastic  and  viscous  matrices;  (ii)  the  first  row  of  the  acoustic  layer  matrix 
and  (iii)  the  lower  half-space  element.  An  example  of  an  assembled  system 
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matrix  is  shown  elsewhere  [1]. 


d.  Spectral  Excitation 


A  point  force  Fq  at  a  layer  interface  may  be  represented  as  an  external 

stress 

E<x,y)  =  Fgetx-Xfl)  Sty-yg)  (2.4) 

whose  'spectral'  form  is  obtained  via  equation  (2.2)  as 


j  E(.a.0)  =  Foexp(-iaXo-ijSyQ)  (2.5) 

For  the  case  of  a  point  source,  of  free-field  pressure  pgexp(  ikRg) /Rq  . 
located  at  distance  H  below  the  lowest  interface,  the  z-direction  'spectral' 
excitation  in  the  absence  of  flow  is  stated  without  proof  as 

! 

E^pta.iS)  =  47ripgexp(-iaxg-i;eyg+i’y  I  H  I )  (2.6) 


e.  Free-Wave  Propagation 


In  the  absence  of  external  forces,  the  system  of  homogeneous  equations 
(2.3)  has  a  non-trivial  solution  only  if  detlZ( a, ^8) ]  vanishes.  For  a  given  value 
of  w,  there  will  be  real  and  imaginary  values  of  a  and  /3  at  which  the 
determinant  vanishes.  Real  values  of  (.a.S)  are  the  wavenumbers  at  which  free 
waves  propagate.  Complex  values  of  (.a.0)  describe  the  evanescent  waves 
whose  effect  decreases  exponentially  with  distance.  Here,  the  real  branches 
only  are  of  interest.  Plots  of  real  wavenumber  a  versus  frequency  for  a 
selected  value  of  wavenumber  0  are  called  wavenumber-frequency  or  dispersion 
plots. 


f.  Far-Field  Pressure 


The  Fourier  transform  representation  enables  the  far-field  radiated  sound  in 
the  upper  half-space  to  be  given  explicitly  in  terms  of  the  'spectral'  normal 
displacement  of  the  uppermost  interface  as  [2] 


p(R,0,<p)  =  -pw^tu^i  (a. /8.  h)exp(ikRg) /27rR]((l-aM/k)  2/q]  (2.7) 
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where 


g=[-M.  sin(  0)  cos(<p) +q]/ ( 1-M2)  ,  q 2  =  1 -M 2[l-sin 2  (  Q)  cos^  ("P)  ]  (2.8) 


and  the  stationary  phase  point  is  given  by 


a  =  kt-M+sin(  0)  costyj) /q]/ ( 1 -m2  ) 
0  =  k.  sin(  0)  sin(ip) /q 


(2.  9) 


3.  THE  ELASTIC  LAYER 


a.  General 


Figure  2(a)  shows  a  section  through  an  elastic  layer  with  lower  boundary 
z=0  and  upper  boundary  z=h.  The  surfaces  are  subject  to  prescribed  normal 
and  tangential  spectral'  stresses 


[T(a.0)]  = 

[T2x(a.  0.  h)  .  r^yt  a.  0.  h) ,  722(0!,  0,  h) ,  0.0).  T2y(a,  0,0),  0,0)]^ 

(3.  1) 

which  produce  'spectral'  surface  displacements 


[u(a.0)l  = 

[Ux(a,0,h),Uy(a,0.h).U2(a!,0,h),Ux(a!.0.O).Uy(a,0.O).U2(a!-0.O)]^ 


(3,  2) 


A  matrix  [S(a.0)]  is  required,  which  relates  the  surface  'spectral'  stresses 
and  displacements  by  the  equation 


(S(a!.0)]lu(a!.0)) 
6x6  6x1 


[7(0!. 0)) 
6x1 


(3  3) 


b.  Equations  of  Motion 


The  linear  elastic  equations  of  motion  [4] 


(  X+M) grad[div( u) -  pd^^/dx^ 


(3.  4) 


may  be  reduced  to  three  wave  equations 

V^F+kL^F  =  0  V^G+kj^Q  =  o  V^H+kjSH  =  0  (3.5) 

by  means  of  the  substitutions 

Ujj  =  3F/dx+8G/3y-0  2H/3x0z 

Uy  =  3F/3y-3G/3x-32H/3y3z  (3.6) 

U2  =  3F/3z+32n/3x2  +  32H/3y2 

Substituting  the  Fourier  transform  representations  of  F.G  and  H  into 
equations  (3.5)  and  integrating  the  resulting  equations  gives 

Fta.^.z)  =  Ai  (a,  ;S)  exp(  io'LZ) +A2(a.  >3)  exp(-i7L2) 

Q(.a.0.z)  -  A3(af,  >3)  exp( i'v-pz) +A4(a,  ;S)  exp( -i-yjz)  (3.7) 

H(Q!,/?,z)  =  A5(a.  >3)  exp(  i-yjz) +A6(fl!. /8)  exp(-i7T2) 
where  A]  to  A0  are  (he  unknown  constants  of  integration. 


c.  Spectral  Displacements 

The  spectral'  displacements  are  obtained  by  substituting  the  Fourier 
transform  representations  of  F.G  and  H  into  equations  (3.6).  They  are 

Uy^^a,0.z)  =  iaF+i/SG-ia3H/3z 

Uy(a,jg,z)  =  i>8F-iaG-i/(33H/3z  (3.8) 

Uz(cr,^.z)  =  0F/3z-a2H-/32H 

F.G  and  H  may  be  eliminated  from  equations  (3.8)  by  the  use  of  equations 
(3.  7)  to  give 

Ux  =  iaA-)  (a,  ^)  exp(i7LZ)  +  iaA2( a. /8)  exp( -i^L^) 

+  i/3A3(a, /3)  exp(  i7j  z)  +  i;8A4(  a.  >8)  exp( -i^jz) 

0)  exp(  i>jz)  -a7TA0(a,  0'>  exp(  -ivjz) 


Uw  =  iiSA-|  (a,  ;S)  exp(  i-yiz)  +  ijSAofa. exp(-i7LZ) 

jl 

-  iaAsC  a, /3)  exp(  iYjz)  “  $)  exp(.-\yjz)  (3.9) 

+>97TA5(a.  0)  exp(  i7jz)  -^S^jAeCa,  0)  Qxp(.-\yjz) 

=  iyi_Ai  ( a.  0)  exp(iyi_z)  -  iyi_A2(a,  0)  exp(-iyi_z) 

-  a^A^(a.0)exp(  iyjz)-  a^AgCa.  0)  exp(-iyjz) 

~  0^A^(a.  0)  exp(  iyjz)  -  0^AQ(a.  0)  exp(~iyjz) 

d.  Spectral  Stresses 

it  is  straightforward  to  show  that  the  stresses 
r^x  =  m(  3u2/ax+aux/3z) 

r^y  =  /x(  au2/3y+duy/3z)  (3,10) 

Tz2  =\div(u)+2fidu2/dz 

have  the  following  'spectral'  representations  in  terms  of  the  unknowns  Ai  to  A5  , 

f^x  =  -2Ma7LAiexp(  i7|_z) +2<ice7LA2exp(-i7LZ) 

-  /iiS7TA3exp(  i7TZ)  +  /i0yjA4exp(-iyjZ) 

+i/ia(yj^-a^-0^)  A^expdyjz)  +ifia(yj^-a^-0^)  AgexpC-iyjz) 


r^y  =  -2(i0yi_A^exp(iyiz) +2u,0yi_A29xp(-iyiZ) 

+  /iaf7TA3exp( i7TZ)  -  /i0t7TA4exp(-i7T^)  (3.11) 

HM0(  yj^-a^-0^)  Agexp  ( iyjz)  +iti0(  yj^-a^-0^)  Agexp  ( -177^) 


T22  =  -(  Xk|_2+2M'yL^^ 

-2i/x7T(  ^5®xp^  *2iu.yj(  A5exp(-i7TZ) 
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Equations  (3.9)  and  (3.11).  evaluated  at  z=h  and  z^o,  give  the  matrix 
equations 


[R(a,^)](A(a.;3)l 
6x6  6x1 


lula.jS)] 

6x1 


(3. 12) 


% 


[P<.a.0)  mca.fi)] 
6x6  6x1 


[f(a!.>8)) 

6x1 


from  which  (A(a..d)]  may  be  eliminated  to  give  the  required  relation  between 
surface  'spectral'  stresses  and  displacements,  viz. 


[P(a,>3)lIR(a,/3)r''[u(a.iS)]  =  lTla./3)] 
6x6  6x6  6x1  6x1 


(3.  13) 
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The  elements  of  the  matrices  [P(a.;3)]  and  [R(a.i3)]  are  listed  in  Appendix  A. 


f.  Special  Case 


In  the  special  case  a=0=O,  the  matrix  [R(a.^)]  is  singular,  so  the  matrix 
inversion  is  not  possible.  However,  it  is  not  difficult  to  show  that  the  required 
'spectral'  stress-displacement  relation  is 


I 


w2p/KL 

sin(kLh) 


cos(kLh)  -1 
1  -cos(ki_h) 


Uz(0.  0,  h) 

u^CO,  0,  0) 


T^^CO,  0,  h) 


T^^CO.  0.  0) 


(3.  14) 


which  as  expected  is  independent  of  the  shear  wave  speed,  cj.  The  'spectral' 
displacements  Ux-  Uy  are  identically  zero. 
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4.  THE  VISCOUS  FLUID  LAYER 


a.  General 


Figure  2(a)  shows  a  section  through  a  viscous  fluid  layer  with  lower 
boundary  z=0  and  upper  boundai  •  z=h  The  surfaces  are  subject  to 
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prescribed  normal  and  tangential  'spectral'  stresses  which  produce 

'spectral'  surface  displacements  [uta.^)].  A  matrix  is  required  which 

relates  the  surface  'spectral'  stresses  and  displacements  by  t'^e  equation 


tSv(a.iS)]Iu(a!.;3)]  =  Vria.B'i] 
6x6  6x1  6x1 


b.  Eoiiations 

The  linearized  Navier-StoKes  equations  of  a  viscous  fluid  [4] 

-grade  p) +MV2u+(/i/3)  gradtdiv(u)  ]  =  pSu/at 
ap/at+pdiv(  u)  =  0 
ap/at  =  c^ap/at 

may  be  reduced  to  the  equations 

v2F+(w2/(c^-4iwM/3p))F  =  0 
V^G+Ciwp//il  G  =  0 
V^H+e  iwp/^)  H  =  0 
p  =  iwpF+(4^/3)  V^F 

pw  means  of  the  substitutions 

Ux  =  dF/dx+dG/dy-d^H/dxdz 
Uy  =  aF/ay-aG/ax-a^H/ayaz 
=  aF/az+a2H/ax2+a2H/ay2 

Equations  (4.3)  and  (4.4)  together  with  the  stress-velocity  relations 

r^x  *  /u.(duz/dx+dux^dz) 

T^y  =  jj,(  dii^/dy+duy/dz) 

T2z  -  -p-(  2/4/3)  div(^) +2M3uz/az 


(4.  1) 


(4.  2) 


(4.  3) 


(4.  4) 


(4.  5) 


enable  the  analysis  to  proceed  as  In  Section  3.  to  give  the  matrix  relation 


[P(a,/9)](A(0!./?)1 
6x6  6x1 


[R(a,>3)][A(a,/?)] 
6x6  6x1 


[T(.a.0)] 

6x1 


[u(a./3)] 

6x1 


(4.  6) 


from  which  [A(a.^)]  may  be  eliminated  to  give  the  required  relation  between 
'spectral'  stresses  and  displacements,  viz. 


-iw(P(a,;S)](R(a./9)]"''[u(a.>S)]  =  [Tia.0)] 
6x6  6x6  6x1  6x1 


(4.  7) 


The  elements  of  lP(a./3)]  and  [R(a!.;S)l  are  the  same  as  for  the  elastic  layer, 
except  that 


=  [w^/ ( c2-4iwp/3/i)  ]-a2_j9  2 


s  {  iwp/ 

=  (ipc^/w) -2m/3 


(4.  8) 


The  special  case  of  the  singular  matrix  a=/3~0  is  also  treated  in  the 
same  way  as  in  Section  3.  except  that  =  w2/(c2-4jw/i/3p)  . 


a.  General 


Figure  2(&)  shows  a  section  through  a  fluid  layer  with  lower  boundary  z=0 
and  upper  boundary  z=h.  The  fluid  flows  with  constant  subsonic  velocity  V, 
parallel  to  the  x-axis.  The  surfaces  are  subject  to  prescribed  normal  spectral' 
pressures 


fp(a.)S)]  =  (p(a,>0.h)  .  p(at.>3.0))''' 


which  produce,  normal  to  the  surface,  the  'spectral'  displacements 


[u2(a,;S)]  =  [uz(a. /5.  h)  ,  u^ta,  0)) 


(5.  1) 


(5,  2) 


A  matrix  [F(a[,>S)]  is  required,  which  relates  the  'spectral'  pressures  and 
displacements  by  the  equation 


[F(a.;9)l[u2(a!.>3)l  =  tp(a.;S)] 
2x2  2x1  2x1 


(5.  3) 


The  linearized  convected  form  of  the  wave  equation 


=  ( i/c2)ia/3t+va/ax]2^ 


(5.4) 


may  be  put  into  its  'spectral'  form  by  substituting  the  Fourier  transform 
representation  of  y?  into  equation  (5.4)  and  then  integrating  the  resulting  equation 

to  give 


<pCa.0.z)  =  A]  (a, /3)  exp(  iyz) +A2(a,  a)  exp( -i7z) 


(5.  5) 


where 


y  =  Vlk^d-aM/k) 


The  condition  of  continuity  of  normal  displacement  between  the  fluid  and  a 
fixed  point  on  the  boundary  is 


d<p/dz  =  auz/at+vau^/sx 


(5.  6) 


from  which,  using  equation  (5.5).  the  displacement  can  be  shown  to  have  the 
following  'spectral'  representation 


U2(a./3.z)  =  -[7/w(  1-aM/k)  JA]  (a. /3)  exp(  i7z) 
+C7/w(  1  -aM/k)  ]A2(  a.  /3)  exp(  -i7z) 


(5.  7) 


It  is  straightforward  to  show  that  the  pressure 
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p  =  -p(  3«p/3t+\/aip/6x) 


(5.  8) 


has  the  following  'spectral'  representation  in  terms  of  the  unknowns  A]  and  A2 


pia.P.2.)  =  iwp(  l-aM/k)  A-)  ( a.  ;8)  exp(  i7z) 

+iwp(  1-aM/k)  A2(at,  0)  exp(-i7z)  (5.  9) 


e.  Matrix  Equation 


Equations  (5.7)  and  (5.9),  evaluated  at  z=h  and  z=0.  give  the  matrix 
equations 


(R(a.)8)](A(a,;S)l 
2x2  2x1 


(Uz(a,>8)l 

2x1 


CP(a!./3)JtA(Qr./8)l 
2x2  2x1 


[p(a.>S)) 

2x1 


(5.  10) 


from  which  (A(a.>3)]  may  be  eliminated  to  give  the  required  relation  between 
surface  'spectral'  pressure  and  normal  displacement,  viz. 


-w2p(  1-aM/k)  2 

" 

cos('yh)  -1 

Uz(a,  0.  h) 

—  - 

p(a.  0 .  h) 

7sin(7h) 

1  -cos(7h) 

.  . 

0.  0) 

p(  O'.  0) 

6.  THE  ACOUSTIC  INFINITE  HALF-SPACES 

a.  Upper  Half-Space 


Kigure  2(c)  shows  a  section  through  the  upper  acoustic  half-space  with 
boundary  z=0.  The  fluid  flows  with  constant  subsonic  velocity  V.  parallel  to 
the  x-axis.  A  relation  is  required  between  the  prescribed  spectral'  pressure. 
pi(X.0.O).  and  the  normal  displacement.  Ujta./S.O). 

The  mathematical  analysis  is  similar  to  that  used  in  Section  5.  except  that  A2 
is  identically  zero  because  outgoing  waves  only  are  allowed.  The  boundary 
condition  equation  (5.6)  and  the  pressure-potential  relation,  equation  (5.8). 


give  the  required  relation  as 


[Fu(a.;3)](Lrz(«‘/^>l  =  (p(cr.;9)] 
1x1  1x1  1x1 


(6.  1) 


where 


Pyta.;?)  =  -ipw2(  i-crM/k)  2/y 


b.  Lower  Half-Soace 


Figure  2(0)  shows  a  section  through  the  lower  acoustic  half-space  with 
boundary  z=0.  The  fluid  flows  with  constant  subsonic  velocity  V,  parallel  to 
the  x-axis.  In  this  case.  A-j  is  identically  zero  because  outgoing  waves  only 
are  allowed.  Thus  the  required  relation  is 


(FL(a,;S)]IU2(«-^>J  =  Ip<a.;S)] 
1x1  1x1  1x1 


(6.  2) 


where 


=  +ipw2(i-aM/k)  2/7 


7.  FORTRAN  COMPUTER  PROGRAMS 


Fortran  computer  programs  have  been  written  for  a  POP- 11/34A  computer 
running  under  the  RT-11  operating  system.  The  programs  calculate  the  far-field 
sound  radiatron  due  to  point-source  or  point-force  excitation,  and  also  the  real 
branches  of  the  wavenumber  a  versus  frequency  plots  for  a  selected  value  of  0. 
A  wavenumber  is  found  by  stepping  through  a  range  of  a-values  until  a  sign 
change  occurs  in  det(Z(a.>9)]:  it  is  refined  to  a  selected  accuracy  by  repeated 
interval  halving. 

Programs  have  been  written,  using  complex  arithmetic,  in  both  singie 
precision  (32  bit  accuracy)  and  double  precision  (64  bit  accuracy).  The  double 
precision  version  was  necessary  because  the  singie  precision  arithmetic  gave 
numerical  problems  of  ill-conditioning  at  the  low  frequencies  which  were  of 
special  interest  in  the  wavenumber  versus  frequency  plots.  Because  double 
precision  complex  arithmetic  is  not  a  standard  feature  of  RT-11  Fortran,  it  was 
necessary  to  write  a  series  of  subroutines  that  simulated  double  precision 
complex  arithmetic  by  operating  on  double  precision  real  arrays  of  leading 
dimension  two.  Details  of  these  subroutines  together  with  some  equivalent 
assembler  level  routines  are  to  be  published  elsewhere  [51. 


The  material  and  geometric  constants,  in  SI  units,  used  in  the 
computations  are  as  follows: 

I  steel  plate:  X=10.44E10.  m=7.56E10.  p=7700. 0.  h=0.  01 

water:  p=1000.0.  c=1500.0 

Damping  was  included  in  the  sound  radiation  computations  by  setting 
X=X ( l-i7j and  where  the  hysteretic  loss  factors  were  chosen  as 

0.  02. 

I 

b.  Wavenumber-Freouency:  Single  Layer 


Figures  3  and  4  are  wavenumber  versus  frequency  plots  of  a  steel  plate 
I  loaded  on  one  side  only  with  water  which  moves  in  the  positive  x-direction  with 

:  velocities  7.  5  and  22.  5m/s.  respectively.  For  comparison  purposes. 

wavenumber  versus  frequency  plots  without  flow  are  also  given.  In  Figure  3. 
there  are  no  numerical  results  below  0.25Hz.  due  to  ill-conditioning  despite 
double  precision  arithmetic  being  used  for  the  calculations. 

j  The  plots  with  flow  show  the  ch.aracteristics  of  instabilities,  viz.  .  infinite 

!  group  velocity  (dflt/dw)"^  in  Figure  4  and  a  multivalued  wavenumber  branch  in 

Figure  3  The  plots  are  in  agreement  with  those  obtained  by  Atkins  12)  who  has 
used  thin-plate'  theory  in  his  analysis  which  shows  that  the  instabilities  are  of 
the  type  convective  (growth  downstream)  and  absolute  (growth  at  a  point), 
respectively. 

I 


c.  Wavenumber-Freouencv:  Three  Lavers 


Figure  5  contains  wavenumber  versus  frequency  plots  of  a  system 
comprising  a  5cm  layer  of  water  which  is  sandwiched  between  two  steel  plates. 
Omitted  from  the  plots  is  a  wavenumber  branch  that  is  close  to  a  fluid  wave  of 
wavenumber  K.  The  plot  with  a  flow  speed  of  22.  5m/s  contains  a  multivalued 
wavenumber  branch,  but  further  work  using  a  root-locus  technique  (2)  would  be 
necessary  to  confirm  the  presence  of  the  instability  and  to  determine  its  nature. 


d.  Sound  Radiation 


Figures  6(a)  and  6(b)  show  the  far-field  sound  radiation  due  to  point  force 
excitation  of  the  lower  facing  of  a  system  composed  of  a  5cm  layer  of  water 
sandwiched  between  two  1cm  steel  plates;  the  upper  and  lower  half-spaces  are 
water  and  vacuum  respectively.  The  sound  levels  are  expressed  in  dB  reference 
1  micropascal  at  Im  for  IkN  rms  force. 


Figure  6(a)  shows  the  sound  level  variation  with  frequency  at  9=0°’:  the  first 
resonant  peaK  is  caused  by  a  coupled  plate-fluid  resonance,  while  subsequent 
peaks  are  close  to  the  frequencies  at  which  'half-wavelength'  transmission  is 
possible  by  plane  waves  at  'normal'  incidence.  Figure  6(b)  shows  the  angular 
variation  of  the  radiated  sound  at  a  frequency  of  40KHz;  the  broad  lobes  centred 
at  9=62“  are  the  well-known  coincidence  lobes  (31  of  a  1cm  plate:  the  small 
peaks  at  9=17“  are  caused  by  longitudinal  waves  in  the  plates;  other  peaks  are 
caused  by  coupled  layer/plate  resonances  corresponding  to  oblique  incidence  of 
plane  waves  -  their  physical  interpretation  via  wavenumber-frequency  plots  is  to 
be  discussed  elsewhere. 
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FIO.  5  WAVENUMBER  ALPHA  VERSUS  FREQUENCY.  BETA-O.O, 
TWO  ICM  STEEL  PLATES  SEPARATED  BY  5CM  WATER. 


The  Elastic  Layer  Matrices 

-2nay\Bxp(.  i'VLh) 

exp( iyjh) 

-2iLBy\_axp(.\y\h) 

/i«7T®3<p(  i^jh) 

(-XkL2-2/4'yL2]GXp( 

0 

-2\^.yJ(.a^■^0^)expi'\yJh) 

Rows  4,  5  and  6  are  obtained  by 


iPia.B)]  and  [R(a.>3)l 
6x6  6x6 

Pi2=  2iLcey\_Qxp(.-\y\J\'i 
Pi4=  /i^Syjexpt-iyjh) 

Pl5=  ®’'P^“''^T^^ 

P2  2=  2M/87Lexp(-i7Lh) 

P2  4=  -AeryT®xpC-i7T^^ 

P26=  ®’‘P^“'TT^') 

P3  2=  t-Xki_2-2M7L2]exp(-i7Lh) 
P34=  0 

P35=  2iM77(a2+j6^)exp(-i7Th) 

setting  h^O  in  rows  1.  2  and  3. 


iaexp(  iyi_h) 

Rj  2  = 

iaexp(-i7Lh) 

i/7exp(  iy^h) 

Ri4  = 

i/8exp(-i7T*^^ 

ayjexpt  irT*^^ 

Ri6  = 

-oryjexpt-iYT^'^ 

i/9exp<  i7|_h) 

^22“ 

i/3exp(-i7Lh) 

-iaexpt  i7T^^ 

^24“ 

-iaexp(  -i7T*^^ 

ByjoxpC  iyjh) 

1^26“ 

-;37t®*P(“'7t^') 

i7L®xp(i7L*’> 

^3  2" 

-i7L®*p(“i7Ll^) 

0 

R3  4  = 

0 

[-a^-B^]exp(  iyjh) 

*^3  S“ 

[_<3J  2  2  ]Qxp  ( -  i7T*l ) 

Rows  4,  5  and  6  are  obtained  by  setting  h=0  in  rows  1.  2  and  3. 
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